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The center of mass of the Earth is the natural and unambiguous origin of a geo- 
centric satellite dynamical system. A geocentric reference frame assumes that the 
origin of its coordinate axes is at the geocenter , in which all relevant observations 
and results can be referred and in which geodynamic theories or models for the dy- 
namic behavior of Earth can be formulated . In practice, however , a Ainernaticaliy 
obtained terrestrial reference frame may assume an origin other than the geocen- 
ter . A fast and accurate method of determining origin ofTset from the geocenter is 
highly desirabie. Global Positioning System (GPS) measurements , because of their 
abundance and broad distribution , provide a powerful tool to obtain this origin 
offset in a short period of time. Two effective strategies have been devised. Data 
from the first Central and South America (Casa Uno) global GPS experiment have 
been studied to demonstrate the ability of recovering the geocenter location with 
present-day GPS satellites and receivers. 

I. Introduction 

Reference frames are established in order to represent 
positions and motions of objects with respect either to the 
Earth (terrestrial frames) or to celestial bodies in space 
(celestial inertial frames). One of the geophysical require- 
ments of a reference frame is that other geophysical mea- 
surements can be related to it; for example, the reference 
frame used for expressing the Earth’s gravity field as a 
spherical harmonic expansion adopts the center of mass 
of the Earth as origin. This frame must be related to the 
adopted terrestrial frame as well as to the inertial frame in 
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which the satellite orbits are calculated. The geocenter is 
at a focus of the orbital ellipse of a geocentric satellite dy- 
namical system. It is therefore directly accessible through 
dynamical methods. But in practice the origin assumed by 
a kinematically obtained terrestrial reference frame can be 
at some location other than at the Earth’s center of mass. 
Any time-dependent ofTset in a geocentric terrestrial frame 
origin from the geocenter can be misinterpreted as plate 
motions. In order to avoid such confusion, it is important 
to determine as accurately as possible the translational ofT- 
set of the adopted terrestrial reference frame origin from 
the geocenter. 

Two effective strategies were devised to determine the 
reference frame origin ofTset from the geocenter using 
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Global Positioning System (GPS) observations. The re- 
sults of covariance analyses performed to investigate the 
accuracy with which the geocenter position can be deter- 
mined with these two strategies have been published in [1]. 
The analyses indicate that the geocenter position can be 
determined to an accuracy of a few centimeters with just 
one day of precise GPS pseudorange and carrier phase data 
collected by a global GPS tracking network. 

Several regional GPS experiments have been carried 
out since 1985, but the tracking ground networks have 
not been extensive enough to provide the global coverage 
needed for accurate geocenter estimation. The Casa Uno 
experiment, however, used a semiglobal network stretch- 
ing over the South Pacific and across the continental 
United States and Europe, thus providing an opportunity 
to demonstrate the capability of GPS data for geocenter 
determination. The analyses of the results using this set 
of data show that the accuracy is strongly limited by a 
nonglobal GPS constellation, the received data quality, 
and the geometrical distribution of the semiglobal ground 
tracking sites. 

II. Two Effective Strategies to Determine 
the Geocenter 

The two strategies devised in [1] to determine the ori- 
gin offset from the geocenter are briefly reviewed here. A 
fiducial network consists of two or more GPS tracking sta- 
tions whose positions have been determined in an Earth- 
fixed coordinate frame to a very high accuracy, usually 
by Very Long Baseline Interferometry (VLBI) or Satellite 
Laser Ranging (SLR). Several GPS receivers at other, less 
accurately known stations also observe the GPS satellites 
along with the fiducial network. The data can be pro- 
cessed simultaneously to adjust GPS satellite states and 
the positions of the nonfiducial sites. The fiducial stations 
established by VLBI provide a self-consistent, Earth-fixed 
coordinate frame; thus the improved GPS satellite orbits 
and the nonfiducial stations can be expressed with respect 
to this coordinate frame to a greater accuracy. The fil- 
ter process is designed so that the baselines between a 
reference site and all other nonfiducial sites are adjusted 
along with GPS orbits and the absolute coordinates of 
the reference site. The correction to the reference site co- 
ordinates infers the adjustment of the geocenter position 
coordinates. This concept has been used in Strategy 1, 
where one or two fiducial baselines are fixed or constrained 
by their a priori uncertainties [1]. The orientation of the 
adopted coordinate frame is defined by the fixed baselines, 
and the absolute scaling can be fixed either by the length 
of these baselines or by the Earth’s gravitational constant 
(GM). Both are known to an accuracy of about one part 


in 10 8 . The absolute scale derived from the fixed baseline 
length allows the coordinate frame thus established to be 
consistent with the VLBI frame of the fiducial baselines. 

In Strategy 2, only the longitude of a reference site is 
held fixed; all other site coordinates are adjusted simul- 
taneously along with GPS satellite states. The absolute 
scale is provided by the Earth’s GM. The geocentric ra- 
dius at a station can be obtained from the adjusted periods 
of GPS orbits and pseudorange measurements. The time 
signature of the measurements defines the latitude. The 
coordinate system thus defined will be an Earth-centered, 
Earth-fixed (ECEF) coordinate frame. 

The covariance study carried out in [1] assumed a full 
constellation of eighteen GPS satellites distributed in six 
orbited planes. A data arc spanning over 34 hours from a 
network of six globally distributed tracking stations was 
also assumed. For Strategy 1, the a priori uncertainty 
for the relative positions of the fiducial sites was assumed 
to be 3 cm. P-code pseudorange and carrier phase data 
noise were assumed to be 5 cm and 0.5 cm respectively 
when integrated over 30 minutes and corrected for iono- 
spheric effects by dual-frequency combination. Carrier 
phase biases were adjusted with a large a priori uncer- 
tainty. The abundance and broad distribution of the GPS 
measurements allow the GPS and station clocks to be 
treated as white-noise processes and adjusted to remove 
their effects on the solutions. Also adjusted were the zenith 
tropospheric delays at all ground sites, which were treated 
as random-walk parameters to model the temporal varia- 
tion. Such models have proven to be effective in remov- 
ing tropospheric errors without heavily depleting the data 
strength [2]. The same network of six tracking sites was 
also used to assess Strategy 2. The estimated quantities 
were the GPS satellite states, the coordinates of all six 
sites except for the fixed longitude of the reference site, 
white-noise clocks, random- walk troposphere parameters, 
and carrier phase biases. 

Data arcs of various lengths were used in the covariance 
analysis. In Strategy 1, at the end of 34 hours the formal 
error in origin offset from geocenter is 4.0 cm (rms of all 
three components). The dominating error is due to the 
assumed error in the fiducial baselines. Any improvement 
in the baseline estimates will therefore directly benefit the 
geocenter determination. As the arc length of the data is 
increased, the error due to data noise is reduced. On the 
other hand, the systematic error due to the fiducial base- 
line errors persists. The situation in Strategy 2, however, 
is different. Here, under the assumptions of the study, 
data noise is the primary error source, which can be re- 
duced by increasing the data arc length. In reality there 
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will be other errors, such as multipath, troposphere mis- 
modeling, and higher order ionospheric effects. The origin 
offset error from the geocenter — as predicted by the co- 
variance analysis using this strategy — -was 2.1 cm at the 
end of the 34 hours [1]. 

III. Methodology for Processing Casa 
Uno Data 

The first Central and South America (Casa Uno) GPS 
experiment was carried out from January 18 to 
February 5, 1988. This experiment was the first civilian 
effort at implementing an extended international (15 na- 
tions) GPS satellite tracking network [3]. Twelve globally 
distributed sites were selected (Fig. 1) to provide improved 
global coverage for the geocenter study. 

The collected data are from seven GPS satellites in the 
constellation, distributed in two orbital planes with a sep- 
aration of approximately 120 deg in the right ascension of 
the ascending nodes. Both strategies are applied to de- 
termine the geocenter location. The analyses present the 
achievable accuracy based upon the available suboptimal 
tracking conditions during the Casa Uno experiment. For 
example, the baselines known to higher accuracies are con- 
centrated only within the continental United States and 
are relatively short compared to the extensive area covered 
by the tracking network; the seven satellites distributed in 
only two orbital planes do not provide a good global cov- 
erage; and the pseudorange data are not of high quality 
due to the antennas and receivers used in this experiment. 

For Strategy 1 the Owens Valley Radio Observatory 
(OVRO)-Haystack baseline is held fixed. The other ground 
station locations, the geocenter location, and the GPS 
satellite states are adjusted with respect to the baseline 
reference point OVRO. Thus the coordinate system has 
the scale and orientation as defined by this fixed base- 
line, and the adjustment to the geocenter location gives 
the offset of the coordinate frame origin from the geocen- 
ter. Table 1 gives the a priori uncertainties adopted for 
Strategy 1. In Strategy 2, the longitude of the reference 
site OVRO is held fixed, and all the other site coordinates 
along with the GPS satellite states are adjusted simulta- 
neously. Table 2 lists the a priori uncertainties adopted in 
Strategy 2, which differ from those in Strategy 1 (Table 1). 

The coordinate systems defined by Strategy 1 and 
Strategy 2 have fundamental differences. In Strategy 1 
the geocenter offset is determined directly while fixing one 
baseline; in other words, the estimated correction to the 
geocenter location represents the coordinate frame offset 
implied by the two ends of the fixed baseline. In Strat- 


egy 2, the coordinate offset is not directly determined. In- 
stead, the corrections to all the tracking site coordinates, 
except for the reference site longitude, are estimated. The 
geocenter offset is inferred from these estimates through 
a constrained seven parameter coordinate transformation. 
The seven parameters are solved by treating filter esti- 
mates of corrections to tracking site coordinates as the 
measurements, whereas the measurement covariance ma- 
trix is represented by the corresponding filter covariance 
plus the a priori covariance of tracking sites. The arbitrar- 
ily fixed reference longitude introduces a small rotation 
R z about the Z axis; this rotation is taken into account 
by solving for R z while the scale factor and the rotations 
about the X and Y axes are kept fixed. The geocenter off- 
set is represented by three translational parameters which 
are assigned a priori uncertainty of 10 m. The Appendix 
gives a full account of the method of seven parameter co- 
ordinate transformation. 

Since the same tracking network was used for both 
strategies, they have similar a priori conditions. For ex- 
ample, in both strategies, the a priori geocentric reference 
frame, in which the nominal station coordinates are ex- 
pressed, is derived from Goddard global VLBI coordinates 
for those same sites and is rotated and translated using 
the results from SLR data for geocentricity. 

IV. Results 

Although the Casa Uno ground network was designed 
to collect a maximum amount of data with the available 
seven GPS satellites, it provides only a semiglobal cover- 
age. In order to compensate for this suboptimal circum- 
stance, the data strength was increased by using a 5-day 
data arc, and the geometrical strength of the network was 
improved by constraining the well-known sites with ap- 
propriate a priori uncertainties. The a priori uncertainties 
assumed for Strategy 1 in this experiment are listed in 
Table 1. In the covariance analysis presented in [1] the 
GPS Block I ROCK 4 model for solar radiation pressure 
was used, where the three parameters G x , G yi and G z 
were considered with 10-percent error. The constant ac- 
celeration in the y axis G y is called the y-bias parameter 
[4], In the Casa Uno multiday analysis a new approach 
was adopted in which two constant solar pressure param- 
eters G y and G xz were estimated along with two tightly 
constrained process noise parameters G x and G z [2]. The 
parameter G xz is a combined effect of ( G x -1- G z ). 

Figure 2 shows the formal rms error in the geocenter 
solution using Strategy 1. The error associated with the 
geocenter offset estimation due to data noise is 31 cm and 
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that due to baseline uncertainty is 7.1 cm, resulting in 
a total formal rms error in the geocenter offset solution 
of 31.8 cm. The effect due to data noise dominates the 
error, showing the poor quality of data; however, the error 
due to the fixed baseline is also significant. These formal 
errors are significantly higher than the 4 cm predicted by 
the earlier covariance study because the ground tracking 
sites were not uniformly distributed, the fixed baseline was 
not long enough to provide a global control, and only one- 
third of the GPS constellation was in place during the 
experiment. The origin offsets estimated with the Casa 
Uno GPS data, with respect to the SLR derived geocenter, 
are —142, —33, and —43 cm along the X, Y, and Z axes 
(Fig. 3). 

Table 2 lists the variation in the a priori uncertain- 
ties from Table 1 as applied to Strategy 2. The longitude 
at OVRO was fixed and all the other stations had large 
a priori sigmas. The same tracking station network was 
used in both cases. The estimated parameters include 
all the ground station coordinates and the GPS satellite 
states. The formal rms error in the geocenter using this 
strategy is found to be 21.5 cm, which is due only to the 
data noise. The estimated origin offsets from the nominal 
geocenter are —96, —97, and —51 cm along the X, Y, and 
Z axes (Fig. 3). 

The origin solutions from the two strategies differ from 
each other along X, Y, and Z axes by 46, 65, and 8 cm re- 
spectively. The rms of these differences is 45.7 cm, which is 
the mean origin offset in each component between results 
from Strategy 1 and Strategy 2. The geocenter estimates 
from Strategy 1 and Strategy 2 agree at the 2-sigma level. 
The apparent disagreement in X and Y components be- 
tween the GPS and SLR solution is puzzling. It may be 
that the baselines between European and U.S. sites are 
known less accurately than the 10 cm assumed here. A 


recent independent (but inconclusive) experiment gave a 
baseline adjustment of up to 50 cm between the station at 
Wettzell and the U.S. sites. 

The anticipated improvement in data quality, better 
GPS constellation, and even distribution of the tracking 
stations should in the future improve the accuracy of geo- 
center estimation to the few centimeters predicted by the 
earlier covariance studies. 

V. Summary and Conclusion 

The Casa Uno 5-day GPS data arc, despite the uneven 
data quality, limited global coverage, and partial satellite 
constellation, can recover the geocenter position to an es- 
timated accuracy of better than 35 cm using either of two 
strategies. The two coordinate frames differ from one an- 
other by three translational parameters and one rotational 
parameter about the Z axis. The transformation parame- 
ters between the two coordinate frames are very sensitive 
to the a priori values and constraints applied to the par- 
ticipating sites. The discrepancy between the results from 
Strategy 1 and Strategy 2 falls within their 2-signia error. 

With future superior receiver data quality, an evenly 
distributed global network including longer fiducial base- 
lines, and an increased number of satellites distributed in 
more orbital planes, the results are expected to improve to 
an accuracy of a few centimeters as indicated by previous 
covariance studies. Inclusion of the Deep Space Network 
(DSN) sites will provide the well determined global base- 
lines that are lacking in the Casa Uno network. Longer 
baseline ties to the well-known sites in North America from 
Europe and the southern hemisphere will provide a much 
stronger network geometry for geocenter determination. 
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Table 1. The a priori assumptions for geocenter study using 
Casa Uno data (one baseline fixed) 


Reference site: 

OVRO 

Other fiducial sites: 

Haystack 

Nonfiducial sites: 

Blackbircli, Canberra, Kokee, 
Samoa, Cocos, Albrook, Mojave, 
Ft. Davis, Wettzell, and Onsala 

GPS constellation: 

GPS 11, GPS 3, GPS 4, 
GPS 6, GPS 8, GPS 10, 
and GPS 9 distributed in 
2 orbital planes 

Cutoff elevation: 

15 deg 

Data span: 

up to 5 days 

Data interval: 

6 min 

Data noise: 

175 cm pseudorange; 
1 cm carrier phase 

GPS epoch state: 

20 km and 20 m/sec (adjusted) 

Geocenter position: 
(each comp.) 

10 m (adjusted) 

Baseline coordinates: 
(relative to OVRO, 
each comp.) 

Haystack 4 cm (fixed); Ft. Davis 
and Mojave — 4 cm (adjusted); 
Wettzell and Onsala — 10 cm 
(adjusted); Others — 1 km each 
comp, (adjusted) 

Solar pressure: 

Gy and G X z (adjusted) 


Table 2. Variation of assumptions from Table 1 for Strategy 2 
(one longitude fixed) 


Reference site: 

OVRO 

Reference site coordinates: 

10 m (latitude) 


0 m (longitude) 


10 m (height) 

Other site coordinates: 

10 m each component 




Fig. 1. The Casa Uno tracking sites for geocenter studies. 



Fig. 2. Geocenter error (1 a) from covariance study for Casa Uno 
5-day arc data using Strategy 1. 



GEOCENTER X GEOCENTER Y GEOCENTER Z 


Fig. 3. Geocenter estimates from Casa Uno data using 
Strategies 1 and 2. 
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Appendix 

Seven Parameter Coordinate Transformation 

The coordinates of a point in a given Cartesian reference frame can be expressed with respect to any other Cartesian 
reference frame as 


x“ 


' AX ' 


X' 

Y 

Z 

Refl 

AY 

AZ 

+ (1 + AL) R 

Y 

Z 


(A-l) 


where AX, AY, and AZ are the three translational parameters, and AL is the scale difference between reference frames 
Refl and Ref2. The rotation matrix R represents the rotation required about the three axes in order to align the two 

reference frames 


R = Yl(R y )Jl(R x )H(R z ) 


where R(/?,) is the matrix representing a right-handed rotation about the t'th axis through an angle Ri. The rotation 
matrix It is obtained by multiplying the three matrices in sequence, i.e., 


R = 


cos R z cos Ry - sin R z sin R x sin Ry sin R z cos Ry + cos R z sin R x sin Ry - cos R x sin R y 
- sin R z cos R x cos R z cos R x sin R x 

cos R z sin R y - sin R z sin R x cos Ry sin R z sin Ry - cos R z sin R x cos Ry cos R x cos Ry _ 


(A-2) 


For small angle rotations and neglecting products of small angles. 


R 


1 R Z -Ry 

-Rz 1 Rx 

Ry —R x 1 


(A-3) 


Substituting Eq. (A-3) into (A-l) and rearranging the terms gives: 

AX - Ry Z Ref2 + Rz X Ref2 + (X Re f2 + Rz^ Rct2 ~ ^Z Ref2 ) AL + (X Ref2 - X Re fl) = 0 

AY + R x Z R ef2 - R-z X Ref 2 + (Y R ef2 ~ Rz^KcVl ~ R X Z Re f2) + ( Y Re f2 “ Y Refl ) = 0 (A-4) 

AZ - R*Y R ef2 + RyX-Ren + (Z R ef2 + ftyX Ref2 - Rx Y Re f2) AL + (Z Re f2 - Z Ref i ) = 0 

This set of equations represents the relationship between two closely oriented, closely scaled orthonormal Cartesian 

reference frames. Every point appearing in both frames will generate three observation equations, where the seven 
transformation parameters (AX , AY, AZ, R x , Ry, Rz, AL) relating the two frames are to be uniquely estimated. A 
unique solution requires two stations with their coordinates known in both reference frames and only one component of 
a third station known in both frames. In practice there will be more than three stations participating and the weighted 
least squares solution gives the solution for the parameters. 

Let Lb represent the observations with corresponding variance-covariance matrix El, so that the weight matrix 
p = (El) - 1 - Each observed quantity is expressed as a function of the parameters in the model 

L a = F(X a ) ( A “ 5 ) 
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Here L a represents the theoretical values of the observed quantities and X a the theoretical values of the parameters. The 
Taylor series linearization gives the observation equation as 

V = Ax T (L 0 — L b ) = Ax T L 

where V is the residuals; A = 5F/3X a |xa=Xo is the matrix of the partials evaluated with respect to the a priori value 
X 0 of the parameters; L = L 0 - L b , where L 0 = F(X 0 ); and x = X a - X Q . The least squares estimate of x is 

x = -(A t PA) _ 1 A t PL (A-6) 

with the corresponding covariance matrix 

E* = (A t PA)- x (A-7) 

If a priori knowledge of the covariance matrix corresponding to the parameters exists, then a constrained least squares 
solution of the parameters is given by 


X = -(P ox + A t PA) _ 1 A t PL (A-8) 

and 

s x = (P ox + A^PA)" 1 (A-9) 

where the parameters are constrained by their a priori weight P ox - 


The partial matrix A of Eq. (A-3) with respect to the seven transformation parameters, when evaluated with an 
a priori estimate of the parameters as zero, would become 


1 0 0 

0 

Z iR e f2 

Y lR e f2 

XiR e f2 

0 1 0 

Z lRef2 

0 

X i Ref 2 

Y lR e f2 

0 0 1 

“ YiR e f2 

X lRgf2 

0 

Z l Ref 2 



1 0 0 

0 

Z n R e f2 

Y n R e f2 

X n Ref2 

0 1 0 

Z n Ref2 

0 

— X n R e f2 

Y n Ref2 

0 0 1 

Y n Ref2 

X n Ref2 

0 

Z n Ref2 


Matrix A will have dimension (3n x 7) for n participating stations. The weight matrix P of the observation is a full 
(3n x 3n) matrix. The vector L is of the form 


L(3n,l) = [(XiR e f2 “ ) , (YlRef2 “ YiR e fi) , (Z 1Ref 2 — ZiR e fi), 


(X n Ref2 - X n Refl), (Y n Ref2 “ Y n R e fl), (Z n R e f2 — Z n R e fi)] 


(A-ll) 


The least squares solution is obtained for the parameter vector 

x (7>1) = [AX, AV, AZ, R x ,Ry,R z , AL] T 


(A-12) 
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